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We investigate the influence of spin-orbit coupling on the Kondo effects in carbon nanotube 
quantum dots, using the numerical renormalization group technique. A sufficiently large spin-orbit 
coupling is shown to destroy the 5(7(4) Kondo effects at zero magnetic field, leaving only two SU(2) 
Kondo effects in the one- and three-electron Coulomb blockade valleys. On applying a finite magnetic 
field, two additional, spin-orbit induced SU(2) Kondo effects arise in the three- and two-electron 
valleys. Using physically realistic model parameters, we calculate the differential conductance over 
a range of gate voltages, temperatures and fields. The results agree well with measurements from 
two different experimental devices in the literature, and explain a number of observations that are 
not described within the standard framework of the SU (4) Anderson impurity model. 



I. INTRODUCTION 

Quantum dots fabricated within carbon nanotubes 
(CNTs) have attracted considerable attention in recent 
years (for reviews see e.g. [1,2]). Such devices, in which 
electrons are trapped within a small, strongly interacting 
region of the CNT by an applied electric field, show re- 
markable electronic transport properties 1,2 and may have 
useful applications in future technology. 3 

CNT-based devices are of particular interest due to 
their doubly-degenerate orbital structure which, com- 
bined with electron spin, generates a wealth of basic 
physics. One such phenomenon is the Kondo effect, 4 re- 
sulting from strong electron interactions within the dot. 
It is observed 5-10 when the device is tuned so that the dot 
has a partially-filled shell of electrons: on lowering the 
temperature, the dot's spin/orbital degrees of freedom 
become strongly coupled 4 to those of the leads, leading 
to a complex many-body ground state with an enhanced 
electronic conductance. 11,12 Understanding this effect is 
particularly important for CNT dots, because the in- 
volvement of both spin and orbital degrees of freedom 
generates an SU (4) Kondo effect that persists to consid- 
erably higher temperatures - up to a few Kelvin (i.e. a 
few tenths of an meV) - than the standard ST/ (2) Kondo 
effect in semiconductor devices. 13-15 Using a range of 
many-body techniques, the theory of the St/ (4) Kondo 
effect is now well established 16-32 and e.g. has been 
shown 31 to be in good agreement with experiments 10 per- 
formed in the absence of a magnetic field. 

Another consequence of the interplay between spin and 
orbital degrees of freedom in CNT devices is spin-orbit 
(SO) coupling. Its existence was beautifully demon- 
strated in experiments 33 on a very strongly correlated 
CNT dot, where it generates a splitting of sequential 
tunneling spectra at finite bias, and kinks in the mag- 
netic field dependence of the Coulomb blockade 'stair- 
case' at zero bias. The strength of the SO coupling was 
measured 33 , and for the device studied found to be of 
order 0.2-0.4 meV, varying somewhat between different 



electron shells. 

Comparing the typical energy scales of the Kondo ef- 
fect and the SO coupling, it is striking that both may arise 
on the energy scale of a few tenths of an meV. Two re- 
lated questions then arise: what effect does SO coupling 
have on the standard SC/(4) Kondo theory, and is the 
Kondo/SO competition seen experimentally? We seek to 
answer these questions in this paper. 

Aspects of spin-orbit coupling in CNTs have recently 
been considered theoretically. The origins of the cou- 
pling itself have been determined from direct microscopic 
calculation, 34,35 showing that while a number of distinct 
spin-orbit interactions arise in principle, the dominant 
contribution is the direct coupling between each elec- 
tron's spin and orbital angular momentum. The effect 
of this coupling on the states of the isolated dot has been 
analysed in detail, 36,37 and the resulting sequential tun- 
neling transport properties (arising when the Kondo scale 
is too small to be seen experimentally 33 ) have been cal- 
culated and compared to experiment. 38 Aspects of the 
competition between spin-orbit coupling and the Kondo 
effect have also been studied 39 via an equation of mo- 
tion decoupling scheme in the U — > oo limit. Within this 
rather crude approximation, 40 the splitting of the SU (4) 
Kondo resonance was examined for finite SO coupling 
and magnetic field, and an orbital Kondo effect found at 
finite field in the two-electron Coulomb blockade valley. 39 

In the present work we consider the two-fold orbitally 
degenerate SC/(4) Anderson impurity model (AIM) in a 
magnetic field, with SO coupling, and study it using the 
numerical renormalization group (NRG) 41,42 backed up 
by simple physical arguments. NRG is ideally suited to 
the problem, being known for similar quantum impurity 
models to provide numerically exact results on the low- 
energy /temperature scales relevant to experiment. The 
model itself also has a strong track-record, a previous 
NRG study 31 of the SU(4) AIM in the absence of a 
magnetic field having shown that the all-important low- 
energy Kondo physics is well reproduced when the bare 
model parameters are fitted to high-energy conductance 
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features such as the Coulomb blockade diamonds. 

The paper is laid out as follows. The model is described 
in sec. II, and the relevant theoretical background and 
NRG technique are discussed in sec. III. The behavior 
of the model in the atomic (lead-uncoupled dot) limit is 
outlined in sec. IV, from which simple arguments arc then 
used to deduce the effect of introducing a finite spin-orbit 
interaction. The main body of the paper is sec. V, where 
we present and discuss the results of NRG calculations. 
We begin by considering the zero-bias conductance, as 
a function of gate voltage, temperature and magnetic 
field. Here we make comparison to the experiments of 
Jarillo-Herrero et al., 8 showing that the orbital splitting 
identified empirically in experiment is readily explained 
by the inclusion of spin-orbit coupling in the model. We 
then turn to a discussion of the finite-bias conductance, 
comparing explicitly to the experiments of Makarowski 
et al. 9 and showing that the asymmetry observed in the 
Kondo peaks at finite bias is also well-described by the 
theory. The paper concludes with a brief summary, and 
a discussion of the applicability of the pure 5C/(4) AIM 
to CNT quantum dots. 



II. MODEL 

The basic model used to describe a CNT quantum dot 
is the S77(4) Anderson impurity model 22,30 ' 31 ' 39 , given in 
conventional notation by 4 

H S u(i) = £k " km + XX ( c lm d m + h - c -) 

k.ra k.m 

The 'flavor' index m takes four discrete values, which 
SU(4:) symmetry reflects physically a combination of de- 
generate spin and orbital degrees of freedom: m = (i,cr), 
with i £ {1,2} denoting clockwise and anticlockwise or- 
bits along the z-direction (major axis) of the CNT, and 
cr £ {Til} f° r the 2-components of electron spin. The fi- 
nal two terms in eqn. (1) represent the isolated dot, with 
orbital energy e and charging energy U = e 2 /C (C is the 
dot capacitance); where n m — d\ n d m , N = J2 m h m is the 
total dot number operator, and = d\ a creates a cr-spin 
electron in orbital i. The first pair of terms describe the 
non-interacting conduction band (lead), and tunnel cou- 
pling between the dot/lead. Each is taken to be spin and 
orbital conserving 22 ' 43 , reflecting physically the fact (see 
e.g. [9] ) that in clean CNT devices the leads are formed 
within the nanotube and so 'carry' the orbital symmetry, 
which is then conserved in the tunneling process. 

The model represents the experimentally relevant sit- 
uation in which the single-particle level spacing of the 
dot exceeds both the tunnel coupling to the leads and 
the intradot interactions. 10 In this case the four-electron 
'shells' of the dot are filled sequentially on sweeping the 



gate voltage V g (oc — e) and, for sufficiently-low tempera- 
tures and source-drain biases, only a single shell need be 
considered at a time. Only the direct Coulomb repulsion 
between dot electrons is moreover included; exchange in- 
teractions are generally weaker 44 and are not necessary 31 
to account for the experimental results of e.g. [10]. 

To the 'standard model' above, we add the coupling to 
the external magnetic field, as well as the key SO cou- 
pling of interest here. For a field B applied parallel to 
the nanotube axis, the Zeeman coupling to the spin and 
orbital degrees of freedom takes the form 38,39 

H B = -Bj> a S?+7o*f) (2) 

i 

where 

§ i = \{n^ ~ (3) 

and 

if = +i(n 1T + n n ) if = -\{h 2] + h 2[ ) (4) 

are the ^-components of the spin and orbital-pseudospin 
operators for orbital i; and where the coupling constants 
are j s = g/iB (with g ~ 2 the electron g-factor) and 
70 = 2^i or b (with /i or ;, the orbital moment). 

The SO interaction obtained from detailed micro- 
scopic calculations is rather complicated 34-37 , but in 
practice only direct coupling between electron spin 
and the ^-component of orbital angular momentum is 
relevant. 34 ' 35 - 38 ' 39 As explained in [38], the SO interac- 
tion is then given by 

4o = -2A^rfsf (5a) 

i 

= A (S' 2 - st) (5b) 

with A parametrising the strength of the SO coupling; 
and which simple form e.g. explains in essence fully the 
sequential tunneling experiments of [33]. 
The full Hamiltonain considered is thus 

H = H su ^ + Hb + -ffso, (6) 

inclusion of Hb + H so lowering the symmetry from SU (4) 
to U(l) x U(l) x U(l) x U(l) (such that H commutes 
with the four charge operators H ki(T + n ia ). 

Finally, to connect to experiment we simply treat U, 
A and the ratio 7o/7s as parameters chosen to fit ex- 
periment (which in practice is quite straightforward and 
unambiguous). The dot level energy e is proportional to 
the experimental gate voltage, of form — e = aV g +(3; but 
in practice it is more convenient to work with a dimen- 
sionless gate voltage 31 

such that differences in V g are proportional to differences 
in A^g (with a proportionality constant determined if de- 
sired by fit to experiment). 
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III. THEORETICAL BACKGROUND 

The dot Green function G icr (oj) <-» G icr (t) = 
—iO(t)({c icr (t) 7 cl a }) is central to understanding trans- 
port through the dot; which is directly related to the 
single-particle spectrum D ia {uo) = — i ImG iff (w) via the 
Mcir-Wingrcen formula 45 , as now briefly summarised. 

One first partitions the conduction band of eqn. (1) 
into two equivalent leads, left (L) and right (R). These 
are taken conventionally 4 to be flat bands of width 2D, 
with density of states p(ui) = po = 1/(2D) for \u)\ < D 
(and with D by far the largest energy scale in the model). 
The tunneling matrix elements to the L, R leads are taken 
for simplicity as Vj, and V R , independent of k. The leads 
are fixed at different chemical potentials, p L and p Rl 
with a bias voltage V s d between them, pl — PR = eKid- 
After transients have subsided, the bias voltage gives rise 
to a steady state current through the dot, J, carried by 
its four conduction channels. An exact expression for J 
follows from the Keldysh formalism: 45 

J = TGon(T L + T R ) / dcu [f L {u) - f R (cu)] V D ia (u) 

(8) 



where 



(r L + r R y 



(9) 



Ti, = Tr\V,s\ 2 po is the hybridization strength of the dot 

to lead v (= L or R), and f v (u>) = [e^~^ + 1)] _1 is 
the Fermi function for the lead with inverse temperature 
(3 = 1/T (&b = 1). It is convenient to define r L +T R = T 
(we later take T to be the 'unit' of energy), such that 
the relative strength of coupling to the L and R leads 
enter through the dimensionless Go (which can be cho- 
sen according to the experimental setup under consid- 
eration). In the perfectly symmetric case of Tl = F R , 
Go = 1 is maximal, while in the extreme asymmetric 
case of T L < T R (say), G - 4T L /T R < 1. 

The key experimental quantity is the differential con- 
ductance, G c (T4d) = dJ/dT4d. In the zero-bias limit, 
gives an exact expression for G c (0) = G° in terms of the 
equilibrium single-particle spectrum: 



with f(u) = (e^ w + 1) 1 , which further reduces to 



(10) 



(11) 



for T = 0. Di a {uj) can be calculated accurately at equi- 
librium using the recent FDM NRG method 46 ' 47 . 

To make connection to experiments at finite source- 
drain bias, note that while eqn. (8) itself remains exact 



for finite V s d, the difficulty lies in calculating Di a (u>) out 
of equilibrium. While recent progress has been made 
in applying NRG to the non-equilibrium single-impurity 
Anderson model (see e.g. [48]) it is currently prohibitive 
to apply these techniques to the model of eqn. (6). As in 
previous work 31,49 we thus make the standard approxi- 
mation of neglecting the y s d-dependence of the dot self- 
energy. Using eqn. (8) this leads to 



G c (V sd ) = 



0/lM df R (co) 



5> 



(12) 



where we have taken a symmetric voltage split 49 between 
the leads, Pl/r = ±eVsd/2. Equations (10)-(12) form 
the basis of our calculations of sec. V. 



A. Friedel sum rule 

Eq. (11) relates exactly the zero-bias conductance at 
T = to the four spectra Di a at the Fermi level, u> = 0. 
These in turn can be obtained exactly in terms of the 
so-called 'excess charges' of the dot in the four distinct 
conduction channels, as now sketched. 

The Green function Gi a (oj) is diagonal in spin and or- 
bital indices, and given by 



G ia (uj) =[uj + i0+ - e ia - T(uj) - 1 



(13) 



where T(uj) = T [tt" 1 In \(oj + D)/(u - D)\ - W(D - \u\)] 
is the w-dependent hybridization function, the effec- 
tive one-electron energy under H (eqn. (6)), and S io .(u;) 
is the dot self-energy. Luttinger's integral theorem, 4,50 



Im 



/X 



dui 



G irT (uj) = 0, 



(14) 



applies separately within each conduction channel 
allowing one to follow the steps of e.g. [4] to derive the 
Friedel sum rule: 4,51 



S il7 = TTU:. 



(15) 



This relates 5i a , the (Fermi level) phase shift of the con- 
duction electrons in the (i, a) channel, to the correspond- 
ing excess charge given by 



ii-, 



J -co I k 



where D-^i a (uj) [D^^iu))] is the (i, a) conduction electron 
spectrum for wavevector k in the presence [absence] of 
the dot. The Fermi level value of the spectrum at T = 
is readily shown to satisfy 4 



7iTA<r(0) = sin 2 (<5 lff ), 



(17) 
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and hence from eqn. (11) we obtain 

^ T =°^Go^sin 2 (7rn inw J. (18) 

i.fT 

The T = zero-bias conductance is thus related to the 
excess charges in the four conduction channels (them- 
selves readily obtained via a thermodynamic NRG calcu- 
lation). In the experimentally relevant limit where D is 
the largest energy scale, these excess charges are more- 
over confined to the dot itself. One can then approximate 
n imp . ia by (n ia ), thereby producing a simple relationship 
between the dot occupancy and its transport properties. 

In the 5J7(4)-symmetric limit (A = = B) con- 
sidered previously, 31 the excess charges are equivalent 
in all four channels and hence n- „.• = n imr ,/4 with 

"imp = E 4!( r«imp;i<T- E Q- ( 18 ) thcn reduces to 31 

^°^osin^). (19) 

B. NRG 

We analyse the model eqn. (6) using the numerical 
renormalization group (NRG). This technique 42 has long 
provided access to numerically-exact results for thermo- 
dynamics and, with the recent identification of its com- 
plete Fock space, 52 an equally systematic and reliable 
route to dynamical properties. 

The basic approach is detailed in e.g. [41]. A loga- 
rithmic discretization of the conduction band states is 
first used to map the Hamiltonian onto a countably in- 
finite one-dimensional chain. The linear chain is diago- 
nalised iteratively, starting from a single site and adding 
the others one by one. The key advantage of logarith- 
mic discretization is that the coupling constants along 
the chain decrease rapidly, and the high-energy states 
of one iteration can be discarded without affecting the 
low-energy states retained in later iterations. As such, 
a fixed number of states can in practice be kept at each 
iteration, rendering the iterative diagonalization of the 
Hamiltonian numerically tractable. 

The information obtained from each iteration allows 
one to build up the thermodynamics and dynamics of the 
model. Eigenstates of a given iteration are used to cal- 
culate thermodynamics at an appropriately-chosen tem- 
perature (sufficiently low that discarded states of earlier 
iterations are unimportant, yet sufficiently high that the 
energy splittings of later iterations are thermally smeared 
out). This effective temperature decreases exponentially 
with the iteration number, and hence allows access to 
thermodynamics on all physical energy scales after only 
a modest number of iterations. 

Dynamics of the model are calculated by means of the 
recent observation 52,53 that the set of all discarded states 
forms a complete basis of the discretized NRG Hamilto- 
nian. By expanding the full density matrix in this basis, 



accurate results for dynamical correlation functions may 
be obtained over a wide range of frequency and temper- 
ature scales. 46,47 Dynamics of the discrete NRG Hamil- 
tonian necessarily arise as a series of isolated poles: in 
order to capture the behavior of the original continuous 
model, one then convolves the discrete spectra with an ap- 
propriate broadening function on a logarithmic scale (see 
e.g. [42]). Potential artefacts of the discretization pro- 
cess are minimised in three standard ways: dot-lead cou- 
plings are prcmultiplicd by the standard factor, 41 ' 54 
the Oliveira 'z-averaging' 55 is used to average discrete 
spectra with different logarithmic discretizations, and the 
Green function is obtained not directly but from its self- 
energy 56 , calculated as a ratio of two correlation func- 
tions where any remaining discretization effects largely 
cancel. 

The calculations in this work have been obtained with 
an NRG discretization parameter A = 3, exploiting the 
full U{1) x U(l) x U(l) x U(l) symmetry of the model. 
We have typically averaged results for five different zs and 
have kept the lowest 2500-4500 states at each iteration. 



IV. ATOMIC LIMIT AND PHYSICAL PICTURE 

Here we show that a physical understanding of the 
interplay between Kondo physics, SO and Zeeman cou- 
plings, follows readily by considering the isolated dot (the 
atomic limit, F = 0).The latter has been considered in 
[38], from which we take results as required; denoting the 
dot states in an obvious notation as | f; — ) = ; — ), 

Hi) = 4|H-)> Hti) = 4r4j-;-)> and so on - 

The energies of the 16 possible isolated dot states fol- 
low directly from eqn. (6) with Vk = 0, and may be clas- 
sified by their total occupation number N — (A) . In the 
absence of SO and Zeeman couplings, A = = B, the 
4!/ [AH (4 — A)!] A-electron states are degenerate, with 
energies £ = Ne + N(N— l)U/2. On switching on A and 
B, the N — and N = 4 states are trivially unaffected, 
while the N = 1 — 3 electron states are split as shown 
schematically in fig. 1. For B = 0, SO coupling splits 
both the A = 1 and N = 3 states into two degenerate 
pairs separated by an energy A. The A = 2 states by 
contrast are split into three groups, of degeneracies 1, 4 
and 1, with relative energies —A, and A respectively. 

Switching on the field further splits the states by both 
spin- and orbital-Zecman effects (fig. 1), and in the A = 1 
sector a singly-degenerate ground state (| t; — )) arises for 
all B > 0. The situation is somewhat more complicated 
in the N — 2 and N — 3 sectors, since in both cases 
competition between SO and Zeeman effects leads to level 
crossings in the ground state. In the N = 2 sector it is the 
orbital Zeeman effect that competes with SO coupling: 
the low-field ground state is | |; J.) as favored by the SO 
interaction eqn. (5b), while at higher fields the ground 
state | IJ.; — ) is favored by the orbital Zeeman interaction 
(eqn. (2)) for the experimentally relevant case 7 Q > 7 S ; 
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FIG. 1: Splittings of the N — 1,2,3 atomic limit states, for 
finite B and A > (relative to their values for B — = A). 

the ground state level crossing occuring at a field 

B Q = — . (20) 

7o 

In the ./V = 3 sector by contrast it is spin Zeeman which 
now competes with SO coupling, producing a level cross- 
ing from the low-field ground state | Tj; J.) to | jj,; f ) at 
a field B s (> B Q ) given by 

B s = -. (21) 

These special values of the field turn out to be central to 
the Kondo physics of the model, as explained below. 

From the atomic limit energies, the ground state 'phase 
diagram' is readily constructed as a function of B and — e 
(oc V g ). Fig. 2 first shows the situation where SO cou- 
pling is absent, A = 0; solid lines marking the boundaries 
between states of different ground state charge. For any 
fixed B, increasing — e generates the familiar Coulomb 
blockade (CB) staircase. When e lies sufficiently in ex- 
cess of the Fermi level, both orbitals are empty; and 
on lowering e through the Fermi level the total number 
of electrons on the dot increases stepwise from zero to 
four. Notice that, for all B, fig. 2 is symmetric under 
reflection about the line — e = |f7 corresponding to the 
midpoint of the N — 2-electron valley (i.e. to replac- 
ing e — > — [e + 3t/]); reflecting for A = the essential 
equivalence of states with N and 4 — N electrons under 
a particle-hole transformation (specifically d\ a d 2 _ a ). 
And at points of degeneracy between N and N + 1 elec- 
trons (fig. 2 solid lines), there is naturally facile zero- 
bias sequential tunneling through the dot (and hence en- 
hanced conductance) when it is connected to the leads. 38 



1 ► B 

FIG. 2: Schematic atomic limit (F = 0) 'phase diagram' in 
the absence of SO coupling, as a function of B and —eocV s , 
showing boundaries between ground states of different charge. 



Degeneracies between states of the same total charge 
do not promote sequential tunneling, but are of course 
vital for the Kondo effects arising from coherent cotun- 
neling processes on coupling to the leads (discussed be- 
low). Such degeneracies arise at zero- field in the N = 1, 
N ~ 2 and N — 3 valleys, where all states of given TV are 
degenerate. For any B > however, there is an immedi- 
ate 'transition' to a singly-degenerate state in each case 
with maximal r z = \ {ni — n^). 

The above situation changes qualitatively on introduc- 
tion of SO coupling; as illustrated in fig. 3, showing the 
generic behavior for A > and j > 7 S (albeit for illus- 
tration using a somewhat smaller ratio "f /"fs = B S /B Q 
than seen experimentally) . This structure is clearly more 
interesting than fig. 2. First, the SO coupling A has a 
dramatic effect at zero field. In both the N — 1 and 
N = 3 valleys the degeneracy of the ground state at 
B = is reduced from 4 to 2 (as in figs. (1,3)); the de- 
gree of freedom associated with this two-fold degeneracy 
being neither a pure spin nor an orbital pseudospin, but 
a mixture of the two. In the N = 2 valley the effect of A 
is even more severe: the 6-fold degenerate ground state 
for A = is replaced by the non-degenerate state | til), 
the other five states again lying 0(A) higher (c/. fig. 1). 

The finite-field ground state level crossings in the 
N = 2, 3 electron sectors (fig. 1), mean as shown in fig. 3 
(dashed lines) that for all e when N — 2 there is a crossing 
from | T; j) to | f J,; — ) at B = B D 39 ; and likewise through- 
out the N = 3 sector, a crossing from | fj.; J.) to | ||; f ) at 
B = B s . Associated with these ground state crossovers 
are naturally kinks in the CB steps seen in fig. 3 (solid 
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FIG. 3: As fig. 2, but for the case of a finite SO coupling A. 
Boundaries between ground states with the same charge are 
shown by dashed lines. 



lines), as discussed in [38] and observed in the sequential 
tunneling experiments of [33] on ultraclean CNT dots. 
By contrast there is no ground state level crossing in the 
N = 1 sector (fig. 1). Hence, as evident in fig. 3, the 
'reflection symmetry' seen in fig. 2 for A = is absent - 
states with N and 4 — N electrons no longer being equiv- 
alent under the particle-hole transformation d\ a <-» d 2 _ a 
(arising because H so (A) — > H so (—A) under such). 

As mentioned above the significance of degeneracies 
between same-charge states is that, on coupling to the 
leads, their associated low energy degrees of freedom can 
be screened by many-body Kondo effects which enhance 
conductance through the dot. 11 ' 12 Two basic Kondo ef- 
fects may in fact arise, 28 SU(4) and SU(2). Only when 
the full model is close to being S'[/(4)-symmetric does the 
former arise (we define 'close to' shortly); the more com- 
mon case, occuring for doubly-degenerate atomic limit 
states where the low-energy effective model maps onto a 
spin-i Kondo model, is the SU(2) Kondo effect. One 
key physical distinction between the two is the low- 
energy /temperature Kondo scale on which they are man- 
ifest, generically denoted Ik- In the strongly correlated 
regime U ^> T, the Kondo scales in the two cases are 4 



SU(2) 



r exp 



T, 



SU(4) 
K 



T exp — — 



nil 
8T 

16T 



(22a) 
(22b) 



modulo prefactors that depend weakly on the bare pa- 



rameters, such that T, 



SU(4) 



K 



SU(2) 



With the above in mind, the essential qualitative 
physics of the model is readily deduced. We start with 
A = 0, and consider specifically the zero-temperature 
limit (the effect of temperature being simply to smear 
out the pristine T = behavior over an energy scale 
0{T)). The appropriate atomic limit picture for A = 
is fig. 2. On coupling to the leads the CB steps, and 
associated zero-bias conductance arising from facile se- 
quential tunneling, are broadened over an energy scale 
O(r). 38 In the N = 1, 2 and 3 electron valleys for B = 0, 
SE/(4) Kondo effects take place 22 ' 25 ' 30 ' 31 (effective low- 
energy SU (4) Kondo, or Coqblin-Schrieffer, models 4 be- 
ing obtained via a Schrieffer- Wolff transformation, on 
retaining the 4!/ [AH (4 — N)l] degenerate dot states in 
the ground state manifold and the cotunneling processes 
that connect them). The Sf(4) Kondo physics in the 
N = 1 — 3 valleys will naturally persist at finite fields 
until 7 S B - 0(Tl U{4) ); while for field strengths m ex- 
cess of the SU(4) Kondo scale, the Kondo effect will be 
destroyed and the conductance will be correspondingly 
low. And for A = 0, the conductance in the N = 1 and 
N = 3 electron valleys as a function of B (or T) will be 
coincident; reflecting the equivalence of associated states 
under the particle-hole transformation discussed above. 

On introducing a finite A, the key quantity in deter- 
mining whether SU (4) or SU (2) Kondo physics prevails 

is the ratio £ = A/T^ U{A) (with T^ U[A) the SU(4) scale 
in the absence of SO coupling) . For £ <C 1 the zero-field 
SU(4) Kondo effects described above are still favorable, 
since the Kondo stabilization energy outweighs the split- 
tings of the atomic limit states in fig. 1. For these small 
As, one thus expects the physics to be essentially un- 
changed from the A = limit. Only when A becomes 
comparable to T^ U{ - A ' will it have a noticeable effect. 

For A 3> by contrast the appropriate starting 

picture is now fig. 3, and SU(A) Kondo effects no longer 
arise. The N = 2 valley will not exhibit any Kondo effect 
at zero-field, since its ground state is singly degenerate. 
The N = 1 and N = 3 valleys for B = will how- 
ever display 5(7(2) Kondo effects involving their doubly- 
degenerate ground states (fig. 3): an effective low-energy 
Kondo model of SU{2) form obviously arises in each case, 
under a Schrieffer- Wolff transformation retaining the ap- 
propriate degenerate pair of dot states indicated in fig. 3. 

For B > 0, the zero-field Kondo effects in the N = 1, 3 
valleys will again 'spill over' into the £?-plane by an 
amount of order their SU(2) Kondo scales (rather less 
than in the SU(A) case, as above). In addition however, 
the level crossings occurring at finite fields in both the 
N = 2 39 and A^ = 3 valleys (fig. 3) means that additional 
SU(2) Kondo effects will now occur at fields B = B a and 
B = B s respectively (the two-fold ground state degen- 
eracy of the free dot states at either field generating an 
SU(2) Kondo model under Schrieffer- Wolff ) . These will 
be discernable as long as they are not subsumed by the 
zero-field SU(2) Kondo effect. This is clearly not an issue 
in the N — 2 valley, no zero-field Kondo effect occurring 
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here anyway for A 3> T^ 11 ^; 



ley it requires (A =) 7 S J5 S » T K 



while in the N 

SU{2 



3 val- 

, readily seen to be 



satisfied since T, 



SU(2) 



< T, 



SU(i) 



K 



In the N — 1 valley 



by contrast, the absence of a ground state level cross- 
ing at finite field (fig. 3) means that the SU(2) Kondo 
effect arising here at zero-field will simply be steadily de- 
stroyed with increasing B, dying out on a scale of order 
7 S B - 0{T^ U{2] ). And since the N = 1 and N = 3 
electron valleys in particular exhibit distinct behavior as 
a function of field, then, as for the atomic limit states 
themselves, the N <-> 4 — N symmetry of the conduc- 
tance as a function of — e oc V g is again absent for A ^ 0. 



V. RESULTS 

The above considerations are purely qualitative, and 
we have analysed the model in detail via NRG, over a 
large parameter space. Here we present a selection of 
results, focusing in particular on parameter regimes ap- 
plicable to experiment. Specifically, we make comparison 
to two experimental works: Makarovski et al. 9 (denoted 
'M') and Jarillo-Herrero et al. s ('JH'). The former device 
is somewhat more strongly correlated than the latter (al- 
though in both cases the ratio U/T is sufficiently large to 
generate non-trivial Kondo behavior), so we can compare 
theory to experiment in two distinct physical regimes. 
Details of how the model parameters are chosen will be 
given at appropriate points in the following. It will also 
be convenient to define and use the reduced parameters: 
U = U/T, A = A/r, f = T/T and B = -f s B/T. 



Zero-bias conductance at B = 



Fig. 4 shows the B = zero-bias conductance as a 
function of the dimensionless gate voltage N g = |(1 — ff) 
(sec. II), for a range of temperatures (T) and three differ- 
ent SO coupling strengths. Here we take an interaction 
U = U/T = 20, so that the relative widths of the CB 
peaks for T 3> Tk are in line with the experiments of M 
(cf. the discussion in [31]). 

Fig. 4(a) shows the A = case, considered in [31]. The 
T = conductance (dotted line) evolves stepwise with 
iV g , and follows eqn. (19) as a function of n imp ; with a 
stepwise increase in n imp itself as the dimensionless gate 
voltage is increased. 31 The latter is a result of the rela- 
tively large U /T, leading to strong charge quantization on 
the dot except when e is within 0(T) of the atomic limit 
charge-degeneracy points at N g — | , | , |, | . On increas- 
ing T the zero-bias conductance is rapidly eroded towards 
the centers of the iV g = 1,2 and 3 electron valleys, re- 
sulting in the familiar Coulomb-blockade valley structure 
with conductance peaks around the atomic limit charge- 
degeneracy points. The temperature scale over which the 
erosion takes place is of course T^- U< " 4 \ given by eqn. (22) 
and obtained numerically 57 as T^ U(A) ~ 0.02r [0.03r] for 




FIG. 4: (Color online.) Zero-bias conductance Gc/(Goe 2 /h) 
as a function of dimensionless gate voltage iV g , at B = 
and for U = 20. Each panel shows the conductance for the 
range of temperatures indicated, with SO coupling strengths 



(a) A = 0, (b) A = 0.02 ~ T } 



SU(4) 
K 



and (c) A 



0.2»T^ (4) . 



N g = 2 [1]. For T^ U{i) < T < T, the HWHM of the 
Coulomb blockade peaks are of order O(r) (not precisely 
T due 31 to electron interactions); while for T>r the CB 
peaks simply broaden to become of width 0(T) instead. 

Fig. 4(b) shows the effect of a non-zero SO coupling of 
order A ~ T^ U ^ 4 \ such that SO coupling competes with 
the SU(A) Kondo effects in the centers of the Coulomb 
blockade valleys. The T = conductance is slightly 
eroded around iV g ~ 2, but is qualitatively unchanged 
elsewhere. Upon increasing T, it is clear that once 
T > A, the conductance appears essentially identical to 
the A = case fig. 4(a), as one expects physically. 

On now considering A 3> T^ U ^ (but still small rela- 
tive to the non- universal scale T), the situation changes 
to that of fig. 4(c). Here the Kondo effect in the two- 
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electron Coulomb-blockade valley at N g ~ 2 is destroyed 
at T — as expected (sec. IV), and the conductance 
remains rather low for all temperatures shown. For the 
N g = 1 and 3 electron CB valleys by contrast, the T = 
conductance is still G°/G a ~ 2e 2 /h (as for the SU(4) 
symmetric limit A = 0, which follows in that case from 



eqn. (19) with n- u 



1 (3) for the center of the Ng = 1 



(3) electron CB valley) . This behavior for large A is now 
however symptomatic of the SU(2) Kondo effect aris- 
ing in that readily understood using the general 
result eqn. (18) for G°, considering the N g — 1 case ex- 
plicitly. Recall (fig. 3) that the free dot ground state for 
N g = 1 is the degenerate pair | f ; — ) and | — ; j) (which 
generate SU(2) Kondo under Schrieffer- Wolff on cotun- 
neling to the leads, sec. IV), for which (nif) = \ = (^24) 



and (fiii) = = {n.2t); an d six 



!mp;wj 



(ni a ) as noted 



in sec. Ill A, eqn. (18) gives directly G°/G ~ 2e 2 /h. 
Note however that although the T = conductance for 
iVg = 1 or 3 barely discriminates between SU(4:) (A = 0, 

fig. 4(a)) and SU(2) (A > T^ U[A \ fig. 4(c)), the erosion 
of conductance with temperature occurs more rapidly in 
the latter case - occurring naturally on the SU (2) Kondo 



scale T, 



SU(2) 



(with T^ U(2) /Y 



8 x 1(T 3 in fig. 4(c)). 



The full T-dependence of the conductance in the cen- 
ters of the CB valleys is shown in fig. 5, G® vs T (on a log- 
scale) for a range of A values, and for N g = 2 (fig. 5(a)) 
and N g — 1 (fig. 5(b), N g = 1 and 3 being equivalent by 
symmetry for B = 0, see e.g. fig. 1). Results for A = 
have been considered in [31]: in each valley the strong 
conductance enhancement due to coherent SU (4) Kondo 
transport is evident in the 'Kondo plateau' for tempera- 
tures T < Tk = T^ U ^ (although note that the univer- 
sal scaling forms of G°(T)/G°(T= 0) in the two valleys 
differ quantitatively, 31 reflecting the two distinct SU(A) 
Kondo effects that arise therein). On a temperature scale 
T ~ U/2, a conductance shoulder is also evident, corre- 
sponding to incoherent sequential tunneling transport. 

In the N g = 2-electron valley (fig. 5(a)) the Kondo 
plateau is progressively destroyed on increasing the SO 
coupling A, as all but the lowest (non-degenerate) atomic 
limit states become projected out of the low-energy man- 
ifold (see fig. 1). For A sufficiently large compared to 

T^- U , a peak is seen to emerge in the conductance at 
T ~ 0(A) (in practice T ~ A/2) - naturally so, this be- 
ing the energy gap to higher SO-split states (fig. 1), and 
mixing in of which enhances the conductance. And in all 
cases shown, the high-temperature (T 3> A) behavior is 
entirely coincident regardless of A. 

As expected from the discussion above, the In- 
dependence of the conductance upon increasing the SO 
coupling A in the center of the one-electron valley 
(fig. 5(b)), shows clearly a crossover from the SU(A) be- 
havior arising for A = 0, to the SU(2) behavior aris- 
ing asymptotically for A 3> T^ 11 ^. This limiting form 
occurs in practice for A > 0.2 (i.e. A/T-^ 17 ^ > 7), 
such that a further increase in A naturally leaves the 
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FIG. 5: (Color online.) Zero bias conductance Gc/(Goe 2 //i) 
vs T for various SO couplings A as indicated, at B = with 
U — 20. For (a) the center of the iV g = 2-electron CB valley; 
(b) the center of the iVg = 1 valley. For B = 0, the behavior 
for AT g = 3 is identical to (b) by symmetry (see fig. 1). 



T-dependence of G° unchanged, as seen in the figure. 



B. Zero-bias conductance at finite B 

We turn now to finite field, fixing the SO coupling A 
and considering G° c = G®(B,N g ) as a function of the 
dimensionless gate voltage N g = |(1 — |^) and field B = 
j s B/T. We first consider T = 0, in terms of which the 
finite-T behavior is readily understood. 

The A = conductance as a function of field 22,26,30 is 
shown in fig. 6(a). At B = it has the stepwise form 
seen in fig. 4, with a maximal conductance plateau of 
G°/Go = 4e 2 //i in the center of the two-electron valley 
(N e = 2) and plateaux of G°/Go — 2e 2 /h in the centers 
of the one- and three-electron valleys. For B ^ 0, the as- 
sociated SU(A) Kondo effects are progressively destroyed: 
four distinct CB peaks instead emerge, centered along the 



G° c /(G e 2 /h) 
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FIG. 6: (Color online.) T = zero-bias conductance maps, 
Gc/(Goe 2 /h) as a function of field B — ■y B B/T and dimen- 
sionless gate voltage N g , for (a) A = and (b) A = 0.2. Bare 
parameters here are U = 20 and Jo/js = 7. 



lines of atomic limit charge degeneracy (cf. fig. 2). The 
SU(4) Kondo behavior at B — persists over a finite 



K 



B-range, but is eventually destroyed for ^ S B 3> Tj 

On introducing a finite SO coupling A = 0.2, the pic- 
ture changes to that of fig. 6(b) (cf. fig. 3). The Coulomb 
blockade lines now show the expected kinks at the fields 
B = B Q and B = B Sl while for sufficiently high B ^> B s , 
SO coupling is of course negligible and the behavior ap- 




B 



FIG. 7: (Color online.) Slices through the N g = 1, jV g = 2 
and iV g = 3 valleys of fig. 6(b), showing the evolution of the 
conductance G° with temperature T as indicated. 



proaches that of fig. 6(a). 

Of primary interest here is the effect of SO coupling on 
the Kondo physics. As seen earlier, at B = the SO cou- 
pling destroys SU(4) Kondo in the two-electron valley, 
eliminating the conductance maximum in this region. It 
also reduces the one- and three-electron zero-field SU (4) 
Kondo effects to SU(2) (the latter apparent in fig. 6(b) 
from a clear reduction in the field strength required to de- 
stroy the Kondo peak). Although SO coupling thus has 
a destructive influence on the Kondo effects for B = 0, 
it leads as discussed in sec. IV to two finite-field SU(2) 
Kondo effects (fig. 6(b)) when level crossings occur in the 
atomic limit ground states. Both lead to a significant en- 
hancement of the T — conductance, of order 2e 2 /h and 

SU(2) 

extending over field ranges 7 s i? ~ T K , resulting in 
marked differences between fig. 6(a) and fig. 6(b). 

The effect of temperature on conductance maps is best 
seen by taking slices through the centers of the CB valleys 
in the (N g , B) plane. The center of the two-electron val- 
ley is (by symmetry) e = —3U/2 for all B, while for the 
one- and three-electron valleys we take the trajectories 



and 



e=-U/2+ \ lo B 



-5U/2 - \ lo B 



(23) 



(24) 



respectively (which approach the centers of the CB val- 
leys in the large- B limit). 
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FIG. 8: (Color online.) As fig. 6(b), but for a reduced coupling 
strength [7 = 6, with A = 1.2 and 7o/7s = 5. 

The results are shown in fig. 7. In the two-electron 
valley the only Kondo effect is the SU(2) Kondo 're- 
vived' at the finite field B = B Q (B a ~ 0.028 here). 39 

SU ( 2) 

The corresponding T K y ' ~ 0.007r, and hence on in- 
creasing T the Kondo effect is in essence destroyed by 
T = T/T = 0.1. Directly analogous comments apply 
to the zero-field Kondo effects arising in the one- and 
three-electron valleys; and, for the latter case, to the ad- 
ditional finite-field SU(2) Kondo effect revived at B = B s 
(B s = 0.2 here). We also note that the clear SO-induced 
asymmetry between the one- and three-electron valleys 
persists even for temperatures T/T ~ 0.1, where the 
finite-field Kondo peak at B s is itself thermally washed 
out: relatively small though it is, the conductance in the 
N g — 3 valley appreciably exceeds that in the -/V g = 1 
valley over a wide B-interval. 



valley brought closer together. 

The results shown in fig. 8 agree well with the experi- 
mental results of JH 8 for the four-electron shell centered 
on V g ~ 3V, 58 as evident from direct comparison with 
fig. 2(b) of JH in the interval 2.5V < V g < 3.5V. And 
more significantly, they provide a natural explanation for 
the observations, as arising from the interplay between 
spin-orbit and Kondo physics. 

The bare parameters employed in fig. 8 are themselves 
consistent with the JH experiment. From fig. 2(b) of JH 
(the C2-D2 line therein) one readily identifies the exper- 
imental B s ~ 3T, and likewise the experimental ratio 
(7o/7s =) B S /B Q ~ 5 (from the C2-D2 and B1-C1 lines). 
Since B s = A/ 7s (eqn. (21)), and 7 s /2 = O.OSSmeVT" 1 
(we take g = 2), the experimental B s gives the SO cou- 
pling constant as A = j s B s ~ 0.35meV - which we note 
is in line with that measured recently in the CNT exper- 
iments of [33] via sequential tunneling spectroscopy at 
finite bias. 

Fitting the atomic limit CB peaks to the JH experi- 
mental data gives t//A ~ 5 and hence U ~ 2meV, in 
good agreement with the heights of the CB diamonds in 
fig. 2 of JH; while the ratio U/T is estimated straight- 
forwardly by comparing the widths of the CB peaks to 
their separation. Finally, the experimental temperature 
T = 0.34K ~ 0.03meV is sufficiently small compared to 
the other scales in the problem, that one can set T = 
with impunity in the NRG calculations. 

As above, we consider the behavior seen in fig. 8 to be 
in striking agreement with fig. 2(b) of JH in the region 
2.5V < V g < 3.5V. We also point out that the experi- 
ment deviates from our calculation above a gate voltage 
V g ~ 3.5V. This arises simply because the levels of the 
adjacent four-electron shell in experiment are brought 
into play by the magnetic field, and at sufficiently high 
B 'interfere' with those arising from the shell considered. 
This is naturally not taken into account in the model 
(although it would be straightforward to incorporate). 



1. Weaker correlations 

Thus far we have focussed on the strongly correlated 
regime where U/T 3> 1. This leads to a pristine sepa- 
ration of energy scales: the CB peaks are separated by 
many times their widths, and the Kondo scales are expo- 
nentially smaller than T. 

On moving to a more moderately correlated regime, 
the energy scales naturally begin to merge, but the essen- 
tial situation remains the same. An example is shown in 
fig. 8, where now U = 6 and A = 1.2 (and 7o/7s has been 
reduced slightly to 5). For U — 6, the resultant 57 zero- 
field SU{2) Kondo scale for example is T^ U(2) ~ 0.5F, 

SU(2) 

and hence T K , A and T all comparable. As seen 
from fig. 8, in comparison to the more strongly corre- 
lated fig. 6 this generates a more pronounced asymmetry 
between the one- and three-electron valleys on increasing 
B, in particular with the CB peaks in the three-electron 



C. Finite bias 

So far we have considered the zero-bias conductance as 
a function of gate voltage, magnetic field and tempera- 
ture. Experimentally there is another 'knob to turn', the 
source-drain bias voltage V s d- As explained in sec. Ill, 
this wc handle approximately using eqn. (12), which re- 
lates the finite-bias conductance to the frequency depen- 
dence of the equilibrium single-particle spectra Di a (w). 

Our main interest here is how finite fields and SO cou- 
pling affect the low-energy Kondo behavior of the con- 
ductance. When A = = B, all four Di a (ui)s share a 
common 59 St/ (4) Kondo resonance in each of the N g = 1, 
2 and 3 CB valleys. On introducing a finite B and A, and 
thus lowering the symmetry to U(l) xU(l)x t/(l) x t/(l), 
each Di a (uj) instead possesses a distinct Kondo reso- 
nance. On application of a field, these resonances shift 
away from the Fermi level w = by different amounts, 
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FIG. 9: (Color online.) Finite-bias differential conductance, as a function of dimensionless gate voltage N g = ^(1 — |j) oc V z 
and source-drain bias eV sd /Y; at T = and for field strengths B of (a) 0, (b) 0.04, (c) 0.2, (d) 0.4, (e) 1.0, (f) 1.4. The other 
bare parameters are U = 20, A — 0.2 and 7o/7s = 5, such that for (b), B = B and for (c), B — B s . 



and at high fields in particular the four resonances are 
sufficiently well separated that the combined spectrum 
J2i a Di a (uj) contains four separate peaks. 22 We also em- 
phasise at this point that, despite occasional naive belief 
to the contrary, the field-induced shifts of the Kondo res- 
onance are not simple 'Zeeman splittings': they have in 
fact a non-linear field dependence brought about by the 
strong interactions on the dot, which can either under- 
estimate or overestimate the Zeeman splitting depending 
on the strength of the field (see e.g. [60] and refs therein 
for a discussion of the SU(2) Anderson model). 

Bearing the above in mind, we consider (fig. 9) finite- 
bias differential conductance maps, G c as a function of 
V s d and gate voltage Ng, with each taken at fixed field. 
Again we start at T — 0, moving to finite- T later when 
comparing to the experiments of M 9 . Taking the limit 
T — » in cqn. (12) shows that G c (V s d) is proportional to 
the average of J2i a Di<j{<jJ = i^eV^d), i-e. the (approxi- 
mate) conductance amounts to a symmetrized combina- 
tion of the total single particle spectrum of the dot. 

The B = = A behavior has been described in [31] 
(fig. 5 therein). Two distinct features arise: the nar- 



row zero-bias ST/ (4) Kondo ridges produced by coherent 
many-body tunneling, and finite bias Coulomb blockade 
diamonds generated by incoherent sequential tunneling. 
The former occur only below Ts of order y^ 1 ^ 4 ) anc } a re 
likewise destroyed by the source-drain bias when eV s d 
becomes of the same order, while the latter are of width 
~ max(r,T) and hence rather more robust. 

Figure 9 shows the behavior for a finite A = 0.2 at 
T = 0, for a range of field strengths. The B — conduc- 
tance is shown in fig. 9(a). The Coulomb blockade dia- 
monds are essentially unchanged from the A = limit 31 
since AcT, and the form of the zero-bias conductance 
is as discussed in relation to fig. 4: on switching on A, 
the one- and three-electron valleys (N g ~ 1 and 3) show 
SU(2) Kondo effects instead of 5Z7(4), while the conduc- 
tance in the two-electron valley at A^g ~ 2 is substantially 
reduced. We now see from fig. 9(a) that the reduction 
of the zero-bias conductance in the two-electron valley 
in fact reflects a splitting of the Kondo resonance in the 
symmetrized spectrum: two narrow conductance peaks 
are seen to arise for A^g ~ 2 when eV^d ~ ±A. 

On slightly increasing the field to B = B a , fig. 9(b), 
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FIG. 10: As fig. 9 but with B = 0.7, f = 0.3, and a colormap 
chosen to be similar to fig. 2 of [9]. The arrows indicate the 
positions of low-energy peaks, see text. Taking F = 0.5meV, 
the eV B d axis extends from —10 to 10 meV, B ~ 3T and 
T ~ 2K, in agreement with fig. 2 of [9]. 



the behavior around eKjd ~ changes. First, the SU (2) 
Kondo resonances in the one- and three-electron valleys 
split (in analogy to the well known behavior of the SU (2) 
AIM in a magnetic field). In the two-electron valley at 
B = B Q by contrast, the orbital SU(2) Kondo effect de- 
scribed earlier arises, and as such the splitting of the 
Kondo resonance here is reduced to zero. 

The next 'special' value of the field is B = B s , shown 
in fig. 9(c). Here the SU{2) spin Kondo effect takes place 
in the three-electron valley, whence the splitting of the 
Kondo resonance seen in fig. 9(b) is reduced to zero. In 
addition, two faint 'shoulders' at a small finite bias can 
just be made out. These are the beginnings of the separa- 
tion of the total spectrum ^2 ia Di CT (w) into four separate 
components at high field (as mentioned above), which we 
discuss in more detail below. 

Finally, for B > B s the atomic limit ground states 
in all valleys are unchanged with increasing B, and the 
low-energy peak splittings in the 14d-dependence of the 
conductance all increase monotonically. Moving from (d) 
to (c) and (f) in fig. 9, we see that once the field becomes 
of order F the Kondo peaks simply merge with the CB di- 
amonds (or the Hubbard satellites in single-particle spec- 
tra terminology) , and the latter themselves begin to split 
from then on. 

At this point we make our first comparison with the 
experiment of M 9 . Keeping the ratio U/T = 20, we take 
r = 0.5meV such that the resulting U = lOmeV is in 
good agreement with the heights of the CB diamonds in 
fig. 2 of M. Our choice of A/T = 0.2 then corresponds to 
A = O.lmeV, which is of the same order of magnitude as 
that measured in another device. 33 And we now consider 
T/T = 0.3 and %B/T = 0.7, to be in line with the ex- 
perimental temperature and field, T — 2K and B = 3T. 
Using a colormap similar to M, and taking Go = 1, we 
obtain fig. 10, which is to be compared with fig. 2 of M. 



Our results are in good agreement with the experiment. 
In particular, we note that the low-energy features iden- 
tified in the experimental paper 9 are reproduced by the 
calculations. In the two-electron valley we see a circu- 
lar region of reduced conductance, while in the one- and 
three-electron valleys we capture two and one low-energy 
peak, respectively (marked by arrows in fig. 10). 

On comparing fig. 10 to the T — results in fig. 9, the 
apparent single peak in the center of the three-electron 
valley is in fact seen to be a consequence of thermal 
broadening: At zero temperature in the three-electron 
valley, a single peak centered at zero-bias occurs only at 
the special fields of B = and B = B s as explained 
above (at any other field, this zero-bias peak is always 
split). While a similar splitting arises also in the one- 
electron valley, it increases monotonically from B = 
and is hence somewhat larger than that of the three- 
electron valley. As a result, the two peaks in the one- 
electron valley remain separate at the particular temper- 
ature f = 0.3 (T ~ 2K) used in fig. 10, while the two 
peaks in the three-electron valley are merged into one. If 
the experiment had been performed at a sufficiently lower 
temperature, we would expect two low-energy peaks in 
both the one- and three-electron valleys. 

To compare further with experiment, we now fix the 
gate voltage to lie at the centers of the one-, two- and 
three-electron valleys (cf. fig. 7) and consider the conduc- 
tance as a function of source-drain bias and field. The 
results are shown in fig. ll(a)-(c) as colormaps, and in 
fig. H(d)-(f) as slices at fixed B = Q, 0.2, 0.4, . . ., 1.6. 
For clarity, the data in fig. ll(d)-(f) have been shifted 
vertically by 0.2(1.6 — 5) to separate the individual lines. 
The figure is to be compared with Fig. 3 of M. 

We see in fig. ll(a)-(c) the evolution of the low-energy 
conductance peaks with increasing field, which again 
agree rather well with experiment. In the one-electron 
valley (fig. 11(a)) the single Kondo peak at B = is seen 
to split into the four spectral features highlighted in M; 
while, as discussed earlier, in the three-electron valley 
(fig. 11(c)) for sufficiently-small B the two lowest-energy 
peaks are merged by thermal broadening. For larger B 
(B > 1 here) the two peaks in the three-electron valley 
do eventually separate in our calculations, as expected on 
physical grounds. This splitting is difficult to see in the 
experiment due to the neighboring four-electron shell be- 
ing brought into play at high field, but should be observ- 
able in a similar device with a larger energy separation 
between shells. 

The destruction of the Kondo effect in the two-electron 
valley with increasing B, as observed in M, is also clearly 
seen (fig. 11(b)), first as a splitting of the B = 
Kondo resonance which then rapidly enlarges to leave 
an almost rectangular-shaped 'hole' in the conductance 
around V s d — (fig. 11(e)). Note also that the data slices 
shown in fig. ll(d-f) are symmetrical about V s d = 0, re- 
flecting the assumption in the calculations of a perfectly 
symmetrical voltage split between the leads (sec. III). 
Their experimental counterparts in fig. 3 of M for the 
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FIG. 11: Field dependence of the low-energy conductance peaks in the centers of the CB valleys, for f = 0.3. (a), (b) and (c) 
show G c /(G e 2 /h) as a function of B and eVad/T, along the centers of the one-, two- and three-electron valleys respectively, 
(d), (e) and (f) show slices of the above data taken at fixed fields B = 0, 0.2, 0.4, . . . , 1.6 (top to bottom). For clarity, the 
data in (d)-(f) have been shifted vertically by 0.2(1.6 — B) to separate the individual lines. 



one- and two-electron valleys (which are not appreciably 
affected by 'overlap' with a higher shell) are somewhat 
asymmetrically disposed about V s d = 0. This can in fact 
be reproduced in calculation by parametrising a small 
degree of asymmetry into the voltage split (without af- 
fecting the essential quality of fig. ll(a-c)), although we 
do not pursue it further here. 



VI. CONCLUSION 

In this paper we have studied the effect of SO coupling 
in carbon nanotube quantum dots, by applying the NRG 
to a modified ST/ (4) Anderson impurity model (AIM). 
Our main focus has been the case in which the SO cou- 
pling is comparable to or exceeds the SU (4) Kondo scale, 
since here the two effects interplay and compete, leading 
to a rich range of physical behavior. The differential 
conductance over a wide parameter space has been cal- 
culated as a function of gate voltage, magnetic field and 
temperature, in order to elucidate the key physics of the 
model. We have moreover shown that the inclusion of 
SO coupling accounts for a number of important experi- 
mental observations in the works of Jarillo-Herrero et al. 8 
and Makarovski et al. 9 , the origin of which stems directly 
from the interplay between SO and Kondo physics. 

To conclude, we comment on the suitability of the 



'pure' SU (4) AIM as a model for carbon nanotube quan- 
tum dots. In [31], experimental data of Makarowski et al. 
at zero-field 10 was found to be in good agreement with 
the SU(A) Anderson model without including SO cou- 
pling, over a wide range of U /T. Given the results of the 
present paper, one naturally asks: why? 

Let us first summarise the experiment. The conduc- 
tance of the experimental device 10 was measured, at sev- 
eral fixed temperatures, as a function of the applied gate 
voltage. The latter was swept over a sufficiently wide 
range that four different electron shells ('Groups I-IV) 
were brought through the Fermi level, one at a time. A 
consequence of varying the gate voltage by this relatively 
large amount was that the tunnel couplings to the dot 
(r) varied from one shell to the next. As a result, the 
Group I data were described by an ST/ (4) model with 
r ~ 0.5meV (and U/T = 20), Groups II and III were 
more consistent with T ~ lmeV (U /T = 10), while Group 
IV had T ~ 2meV (U/T = 5). 

To understand why these data could be described 
by the ST/ (4) model, we note that (a) the experimen- 
tal 'base' temperature was 1.3K ~ O.llmeV and (b) 
the Kondo scales for Groups II-IV are all in excess of 
4K ~ 0.36meV (see fig. 3(b) of [10]). Assuming the SO 
coupling to be comparable to the value O.lmeV consid- 
ered above, Groups II-IV can then be described by a 
pure SU (4) Anderson model (for essentially all temper- 
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atures), since each has a T K appreciably in excess of 
the SO coupling. While the latter is not the case for the 
Group I shell 31 of [10], data were only obtained down to 
a temperature of order A where, according to sec. V A, 
the effects of SO coupling cannot be seen in the zero-field 
behavior alone. Only by examining the behavior of the 
Group I shell in a magnetic field (as considered here), or 
by measuring its zero-field conductance down to a rather 
lower temperature on the order of ~ 0.1K or so, can the 
effects of SO coupling be observed. 
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